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t^. Abstract 

O ■ 

• Physics, as known from our local, around-earth experience, meets 
some of its applicability limits at the time just preceding the period of 

' primeval nucleosynthesis. Attention is focussed here on the effects of 

■ the nucleon size. Radiation-belonging nucleons are found to produce 

q-( an extremely high pressure at kT ~ some tens or hundreds of MeV. 

^ ' Quark deconfinement at higher energies would not change the results. 

1 Introduction 

The standard procedure of Physical Cosmology is to take present-day knowl- 
edge and data and travel backwards in time, applying as well as possible our 
local, laboratory- and observatory-tested Physics. That Physics, as we know 
it today, is able to explain so many of the progressively distant and red- 
shifted data is the best mark we have of its astounding range of validity. The 
remotest time for which we have reliable results is the nucleosynthesis epoch: 
well-established Physics is able to give a fair account of the cosmological 
origin of the lightest elements. We shall here be interested in the period just 
preceding that nucleosynthesis era, which we shall call "pre-nucleosynthesis 
period" (PNS period). 

The end of that period - the beginning of the nucleosynthesis age - must 
correspond to a temperature kT of the order of the deuteron binding-energy, 
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which means a few MeVs and a red-shift zk2x 10 10 . The PNS period could 
also be called the "close-packing period" . The "closely-packed" constituents 
referred to are protons or, more precisely, nucleons. A rough estimate gives 
for their concentration a value around 10 23 cm~ 3 and for their mean free 
path A « (n fo cr) _1 « (10 23 x lCT 26 )" 1 = 1000 cm. The usual approximation 
assuming ideal fluids fails, but the physical assumptions on the large ratio 
between the total volume of the system and the total volume occupied by 
the constituents are valid and current Physics can be expected to hold. 

Protons are remarkably stable (lifetime larger than 1.6 x 10 25 years [0]), 
and neutrons decay into protons. Thus, we can safely suppose that the nucle- 
ons present today in the Universe have been around from the "beginning". 
The values of the critical density and baryon density [see below, equation 
([12]) 1 imply a remnant nucleon density in the range 0.059 < < 0.296 
(nucleon x m~ 3 ) at present time. Now, each nucleon occupies a volume of 
the order 2.2 x 10~ 45 m 3 , which means that at z ~ 10 15 they attain a tightly 
packed state: one nucleon per nucleon volume. This will define for us the 
beginning of the PNS period. The mean free path is then of the order of the 
size of the constituents. In a nutshell: the causally-related Universe has a 
volume Vu ~ 10 81 cm 3 and contains Njy ~ 10 74 nucleons. The total "internal" 
volume of these nucleons is V n ~ 10 35 cm 3 . The Universe had that volume 
when z ~ 10 15 . The assumption of infinite system volume - which underlies 
the thermodynamic limit, as well as the very definition of cross-section - is 
then at least doubtful, and the ideal fluid hypothesis is clearly untenable. 
We shall later refine this crude estimate, but the result will be, not quite 
surprisingly, essentially the same for the remnant protons. The interest of 
the more refined approach rests on its formulas, which can be applied to the 
protons belonging to the radiation bath. 

The PNS period runs consequently between z « 10 10 and z ~ 10 15 . The 
strategy to be followed will be rather circular. Quark deconfinement will be 
ignored to start with and protons will be taken as stable. We shall describe 
them by a potential and find that the pressure related to present-day matter 
tends to an infinite value. That would happen, however, at energies for 
which the notion of potential does not apply and for which deconfinement 
is quite possible. We then reconsider the question from the point of view 
of the radiation-belonging nucleons, and find the same effect at much lower 
energies, for which potentials do have a meaning and there is no possibility 
of deconfinement. 

In a first dealing with such unusual conditions we shall feel justified in 
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taking a naive approach. Instead of facing the intricacies of the high-density 
matter equation of state , we shall content ourselves with reasonable order- 
of-magnitude estimates. The proton will be considered as a hard constituent, 
occupying an irreducible hard-core volume ~ lfermi 3 ~ 10 -39 cm 3 , repre- 
sented by a hard-sphere potential. Potentials have been used from time to 
time in Cosmology, for instance to show how the initial singularity can be 
avoided 0. They should, of course, be carefully handled in relativistic con- 
ditions. We shall be attentive to the energy conditions under which the very 
notion of potential can loose its validity. 

Section is a summary of the Standard Model, actually a commented 
formulary devoted to fixing notation, showing the numbers we use and sum- 
ming up some observation values relevant to our subject. We shall ignore 
non-standard possibilities, as eventual "dark" constituents, and accept usual 
reasonable assumptions, such as the Debye screening which renders electro- 
static effects negligible. Such a review of well-known topics is necessary to 
show how and when we part from the standard procedure. In section |3] a gen- 
eral overview of physical problems appearing in the PNS period is given. In 
particular, we present our assumption that protons keep their identity in en- 
ergies much higher than usually supposed. We then proceed to a discussion of 
the hard-sphere potential and to a naive application to the Friedmann equa- 
tions. The result is that an infinite matter pressure would block the backward 
progress at around z m 10 15 if we take into account only the "remnant" pro- 
tons existing at present time. At those red-shifts the nucleons are relativistic 
and the idea of a hard-sphere potential is unrealistic, but it is easier to ar- 
gue starting from the consideration of the remnant present-day protons. We 
show then (section [5]) that the protons appearing through pair-creation from 
the radiation background produce the same effect at much lower energies. 
Pairs of photons with energy barely enough to produce proton-antiproton 
pairs will create non-relativistic protons and anti-protons, and for these the 
notion of a potential does make sense. Pair creation is a very efficient process: 
the number of created protons is very large already at energies much lower 
than 1 GeV. The pressure blockage, in consequence, takes place at rather 
low red-shifts. The possible meanings of these results are discussed in the 
last section. A briefing on relativistic quantum gases is given in Appendix 
A. 
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2 The standard model 



The large scale evolution of the Universe is described 0, by the two Fried- 
mann equations for the scale parameter a(t): 
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which, once combined, lead to the two equivalent expressions 
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This equation can be alternatively obtained from the vanishing of the covari- 
ant divergence of the source energy-momentum tensor, and reflects simply 
energy conservation. Notation is hopefully obvious: p = e/c 2 is the source 
energy density in mass, p the pressure, A the cosmo logical constant. The 
index "0" will indicate present-day values: the red-shift z, for example, is 
given by 

l + Z = a M- (5) 



a{t) 



The Hubble function 



H(t) 



d(t) d 



\na(t) 



(6) 



a(t) dt 

has the present-day value H = 100 h km sec -1 Mpc~ l = 3.24 x 10~ 18 h sec -1 
(with the parameter h, of the order of unity, encapsulating the uncertainty 
in present-day measurements). 

Given an equation of state in the form p = p(p), equation (f|) can be 
integrated to give 

g3 J eQ e+p(s) 



1 + z 



a(t) 



(7) 



For example, for a pure radiation content the equation of state is p = |e, so 
that 

e z = e Q (1 + zf . 
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The energy density of dust matter, with p = 0, will behave according to 

e z = e (1 + zf . 

Recall that Eq.(^) is a mere consequence of energy conservation. These 
results are independent of the parameters k and A. The relationship between 
z and the Hubble function is easily found. Taking the time derivative of 
Eq.([5|) and comparing with the last expression, one arrives at 4f = — H(t)(l+ 
z), which integrates to 

- f* H(t)dt /_ N 

1 + z = e Jt o y '. (8) 
The critical mass density is 

Pent = = 1-88 x 1(T 26 h 2 kg x m" 3 . (9) 

07TG 

The baryon concentration and mass density are given by 

n b = 11.2 tt b0 h 2 (1 + zf m' 3 (10) 

p b = 1.88 x l0- 26 tt b0 h 2 (1 + zfkg m~ 3 , (11) 
where the parameter f^o = = - fi6a - has observational values in the 

Pcrit 

range 

0.0052 < fi w ^ 2 < 0.026 . (12) 

These values for the critical and baryon density lead to the remnant nucleon 
density range used in the Introduction. In terms of H(t), the equations can 
be written as 

* = 2ffi,-» + % d3) 
d{t) = Hit) a{t) ; f t =-3H(p+^). (15) 



Introducing VL m = — — , = tfft and Vtdt) = equation (O) takes 

pcrit ^ Q a 

the form 

H 2 p kc 2 Ac 2 

-"o Pcrit a n 6rt 
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Notice that Vt m refers to the total amount of content: if baryons and radiation 
are to be considered, then Q m = Q b + fi 7 . The above expression gives on 
present-day values the constraint 

f2 m0 + ^fcO + — + ^70 + ^fco + = 1- (17) 



We have used Qko = — with ). There is a recent evidence 

a Q n 

for a large value of and a small value of f2feo @. Choosing for time and 
length the convenient units 

H l = 3.0857 x lO 17 /^ 1 sec ] 4r = 9 - 25 x lO 25 /^ 1 m , (18) 

the Friedmann equations acquire the simpler forms 

* = (19) 

p crit a[ty 3 

3 „ 9 3 p A 1 k . . 

H = - - H 2 ^— + — . 20 

2 2 c 2 p crit 2 2 a(t) 2 v ' 

Gases at high energies in the presence of pair production are also better 
considered in adapted units, which we introduce here while leaving details to 
Appendix A. First of all, given a particle of mass m, it is convenient to use 

kT , 

as the temperature variable. There are also two lengths of major interest, the 
Compton length and the thermal wavelength. The static Compton length is 
a most natural unit of length: 

Ac = • (22) 
mcr 



For the electron and for the proton, respectively, A e = 3.81 x 10 cm and 
X p = 2.08 x 10~ 14 cm. A natural volume cell for the proton will be A 3 
= 9.0 x 10- 42 cm 3 . 

If (3 — 1/kT is the inverse temperature, (the cube of) the relativistic 
generalization |7| of the thermal wavelength is given by 



A 3 T ([3) = 2 vr^mc 2 6 ( — 2 ) = — XI . (23) 

K 2 (pmc z ) \mc z I r a 2 (1/t) 



-/W / hcY 2 7T 2 e- l ' T 



6 



Here K 2 {x) is the modified Bessel function of second order, whose asymptotic 
behaviors, leading to the nonrelativistic and the ultra-relativistic limits, are 
given in Appendix A. Here we only remark that the non-relativistic limit 
gives the usual expression 

A NR W = \ = h^= ^Xc. (24) 

For instance, a proton at kT w 4 MeV will have A^r « 40 A p . A proton will 
"occupy" a degeneracy- volume A 3 = L42 * 3 /2 — cm 3 , from which every other 
proton will be statistically excluded by Fermi repulsion. The ultra-relativistic 
limit will be 

A UR ((3)=n^ (3hc = tt 2 / 3 ^. (25) 

T 

The pressure p and the mass density p in ([I]) and (0) are those of matter 
and radiation present in the Universe, introduced through their expressions 
for ideal gases. Interactions are only taken into account through reactions 
supposed to take place in restricted conditions. As examples, a weak Thom- 
son scattering lies behind thermal equilibrium before recombination, and pair 
production will be responsible for the existence of a huge number of electrons 
when kT is higher than w 0.5MeV. 

The values A = 0, k = lead to very simple solutions and are helpful 
in providing a qualitative idea of the general picture. They will be used 
as reference cases. We shall later exhibit the expression of H(z) in the pre- 
recombination period, as well as the implicit expression of a(t). Nevertheless, 
in order to get a firmer grip on the relevant contributions and the role of 
each term, it is useful to review the customary discussion on the matter- and 
radiation- dominated ages. 

2.1 Matter— dominated age 

Always in the standard approach, baryons (essentially nucleons) dominate 
the energy content at present time. This domination goes back to the "turn- 
ing point" time given below (equation (f48|) ), when radiation takes over. Pro- 
tons are non-relativistic during all this period. The standard argument runs 
as follows. Matter pressure has the expression p b = n b kT b . It appears, how- 
ever, always in the combination pb + Pb/c 2 = nb[m + kTb/c 2 } = ^[mc 2 + kTb}. 
Thus, pb is negligible in the non-relativistic regime. Putting p = 0, the 
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reference case with A = and k = has for equations 

3 p b 



H = H 



2 pcrit 



The general solution is 



' - t + C. (26) 



H{t) 2 

Let us quickly examine 3 models, two with k = 0, A = and a third case 
with A ^ 0. 

a) Matter— domination: dust Universe 

This unrealistic model supposes matter domination all the time along. It 
takes at the "beginning" H t=0 = oo. The integration constant C vanishes 
and the solution is simply 

H{t) = - . 
v ' 3t 

This means that 

da 2 dt 
a 3 t 

The expressions relating the Hubble function, the expansion parameter, the 
red-shift, and the density follow immediately (we reinsert Hq for conve- 
nience) : 

H 2 



[l + zf. (27) 
a{t) /t\ 2 / 3 



H 2 



- 



a(t ) \t J 
/tn\ 2 / 3 / 2 \ 2/3 

Pcrit Hq 

The age of the Universe can be got from (]28f) , by putting z = 0. One obtains 
t = 2/(3H ) ~ 6.5 x 10 9 years, a rather small number. There is, as said, a 
serious flaw in this exercise-model: it supposes that matter dominates down 
to t ~ 0, which is is far from being the case. Furthermore, the protons cannot, 



- EH n b0 = (i + z y n b0 . (29) 

o 
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of course, be non-relativistic at the high temperatures of the "beginning" and 
matter pressure should be added. Let us see two more realistic cases. 



b) Matter— domination: present time 

Let us go back to the general solution (|26|) and fix the integration constant 
by the present value 

The solution, now more realistic, will be 

Hit) = - TT7~} r. (30) 

W 1 + lH (t -t ) K J 

Integrations of ^ = H[t)dt leads then to 

l + Z = W) = [l + lH (t-t )]^- (31) 

Equation (^) keeps holding. Using that equality, and the value (|9]) of the 
critical density, we find 

p b = 1.878 x KT 26 (1 + zf tt b0 h 2 [kg m~ 3 ] (32) 

Dividing by the proton mass, the number density is 

n b = 11.2 x (1 + zf tt b0 h 2 [m~ 3 ] (33) 

Actually, Qbo = 1 in the reference case we are considering. This gives a few 
nucleons per cubic meter at present time. The age of the Universe is basically 
the same as that for the dust Universe: we look for the time t corresponding 
to z — > oo, and find t — t — 2/(3H ). Equations (|30D and ([31] ) hold from the 
turning point down to present times (provided k = and A = 0). 



c) Matter— domination: k = but A ^ 

Recent evidence for k = and a nonvanishing cosmological constant at 
present time gives to this case the prominent role. 
Let us insert (§) into (p9|), to get 

- 3 f* H(t)dt 

p = Pq e Jt o 
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and then insert this expression into the Friedmann equation 

„2 



H 2 



The time derivative gives 
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- 3 f* H(t)dt Ac 2 

Po e + — 



dH 3 /Ac 2 
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H 



Integration leads then to the expression 
H(t)- 



'Ac 2 



:H/^(t-to) 



-U/^-^o 



3\/^(*-*o) 



(34) 



This expression for H(t) gives H = H Q when £ — > t and tends to (j30|) when 
A — > 0. To have it in terms of more accessible parameters, we may rewrite it 

as 



H(t) = H Jn 



1 



(v^I + l) e 3 ^Vni(t-to) + ( V /H A ~ _ i 



(35) 



To get the relation with z, we notice that 
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2 /4ttG\ 
= 2 — pnfl 



Poll + 2)° + 



Ac 2 



fi 6 (l + 20 3 + n, 



gives 



if 2 



Ac 



r = 2 (ifJ m(1 + z) 

Ac 2 (AtiG\ 
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H 2 
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which together imply 



Alternatively, 



H 



2 _ Ac^ 

3_ 

tt2 Ac' 2 
^0 ~ ~ 



(1 + z) 3 
+ Q\ 



:i + z) 



H — Hq Qa 
1 . 



(36) 
(37) 



10 



It remains to use (|3~4] ) to obtain 



Ac 



2\ 1/3 



1 + z 

\ 6 ) 

which is the same as 
a 



+ 



/ Ac 2 TT 



2/3 



, (38) 



a(t) 



1 + z = (4 fi A ) 1/3 



^ov^t-to) 



(v^a + l) e 3H °v^I(*-*o) + (v^a- 1 



2/3 



• (39) 



The neglected matter pressure is given by 
Pb 



C Pcrit 



9.2 x l(T i4 T b (1 + zf a 



bo ■ 



(40) 



This will be one of our main points. Matter pressure is neglected in usual 
treatments for the reasons given above, which assume an ideal gas. We intend 
to take into account interactions between nucleons, and shall find indications 
of a very abrupt raise of this pressure during the PNS period. 



2.2 Radiation— dominated age 

For photons, the temperature behaves as a frequency so that, by the very 
definition of red-shift, T 7 = T 7 o (1+z). For example, hydrogen recombination 
takes place at T 7 = Tf, « 30007^. This, together with the present value 
T 7 o ~ 2.7 for the thermal background, gives a red-shift (1 + z) ~ 1.1 x 10 3 . 
The mass-equivalent density is therefore 

fK - 2.i x io- 7 r 4 r 2 = i.i x io^ 5 (i + # r 2 = n 70 (i + 2 ) 4 , (4i) 



Pcrit 

where 

fi 7 o = 1.1 x 1(T 5 hT 2 . (42) 

Consider again the reference case k = 0, A = 0. Because e 7 = 3 p 7 and 
p 7 = t-, we have 

H = - 2 H 2 = - 2 -tL- , 

Pcrit 
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leading automatically to 



H 2 = 0,0(1 + zf . (43) 

Solving the equation is only necessary to fix the relation between the time 
parameter and the red-shift. The solution, 

Hit) = — , 
w 2t 



implies 



0,0(1 + z) 2 «o 



^ - fi#V2f . (44) 



Before recombination (that is, for higher z's) there is thermal equilibrium 
between matter and radiation, because electrons are free and the mean free 
path of the photons is very small. An estimate of the energy per photon at 
a certain z can be obtained from fcT 70 = 2.3 x 10~ 10 MeV , which leads to 

A;T 7 = A;T 70 (1 + z) = 2.32 x 1(T 10 (1 + z) MeV. (45) 

For example, an energy of 4 MeV corresponds to z ~ 2 x 10 10 . Thus, the 
thermalized state before recombination makes of T 7 , or its corresponding red- 
shift, the best time parameter. We shall retain for later use the expressions 

P~f 1 P~f — >7 m-8 T>4 L-2 _ ^70 (-> . \4 



7.x 10-^^ = ^(1 + ^ (46) 

Pent d Fcrit d 

^ = 2.31 x 1(T 5 (1 + z) Vl^ . (47) 

At recombination time, ^ ~ 39.5 O - The scale parameter, and conse- 
quently the red-shift, behave quite differently in a matter-dominated Uni- 
verse (|2~q) and in a radiation-dominated one fl44|) . Equation shows that 



radiation becomes dominant at high z's. The turning point, or change of 
regime, takes place when p 7 ~ p b , or 

1 + z ~ 4.3 x 10 4 n w . (48) 

This corresponds to t ~ 7.4 x 1(T 8 Q b0 3/2 /H = 2.2 x 10 10 Q^ /2 sec. When 
there is no thermalization, non-interacting matter pressure is negligible with 
respect to radiation pressure by a factor ^ ~ 
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With our proton-related variables, the temperature during the thermal- 
ized period preceding recombination will be 



r= _I = 2.5xlO- 13 (l + z). (49) 



fcT, 
mc 

In terms of those variables we shall have, for example 



n 7 = 3.8 x KT 39 (1 + ^) 3 A; 3 = 0.24 (j-^J (50) 

and 

n b = 6.3 x 10~ 9 n b0 h 2 (^j . (51) 
2.3 The Friedmann solutions 

Let us now go back to the general equations (|i~9"D and (p0|) . Thermalization 
at z > 10 3 has important formal consequences. Dependence of a single 
temperature makes it more convenient to use the variable z, and much of 
the discussion can be made in terms of energies, by using (|4"5|). Using the 
units given in Eq . (|i"8|) , adding matter (|2~9|) and radiation (f4T|) densities and 
extracting from ([T7|) the value of Eq. (^) becomes 

H\z) = fi 70 (l + ^) 4 - ^ (1 + A 2 + £ + (1 + ^) 3 (1 + 4 - £ " ^70) • (52) 

This will be modified when interactions between nucleons are taken into 
account (see equation fl6"4"D below). In the reference case, 

H 2 = 0,0(1 + zf + (1 - Q j0 ) (1 + ^) 3 , (53) 



which reduces to (|27|) when the last term dominates the right-hand side, and 
to ( P^ ) when the first term dominates. 

All this can be see from another point of view. In fact, we have (1 + 
z)H{t) = — 4|. We can consequently introduce the function f{z) = , 
with /(0) = 1 and a(z) = j^. In the units (|T8|) , equation (|20|) becomes 



;i + ,)| = 3/ + 3-^ + A (1 + ,)2_ A . (54 ) 
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Use now x = 1+z, with f(x = 1) = 1. As radiation pressure largely dominates 
at all time, the only contribution to p will come from ((46|). We shall take 
3 ~2~ — = Q^qx 4 , with fL as given in (H3), and rewrite the equation as 

x ^ = 3f + n^ x 4 + -^x 2 -A. (55) 

The solution is 

k A k A 

f (x) = tt l0 x A ~ — x 2 + - + x 3 (l + — - — - fi 70 ) , (56) 



which is the same as fl52|). This point of view will be of interest when the 
pressure of matter is added. This is due to the fact that pressure appears 
only in the equation (|TJ]) for the derivative of H and not in the expression 
(HI) for H 2 . 

Things are not that simple for the expansion parameter, or for the relation 
between z and t. Noting that x = a /a, we have 



da , A „. A 

a — = W iljoo-o ~ ka 2 + — a 4 + ag(l - ll l0 + k - —) a, 

whose solution is given by 

, , f a{t) ydy f . 

t = t + . (57) 

/a(to) ln i0 a 4 - ky 2 + Ay 4 /3 + y (l - fi 70 + ka^ 2 - A/3) a 3 



Most aspects for z > 10 3 can be discussed in terms of the red-shift, and we 
shall have little use for the time variable. 



3 Microphysics at the PNS period 

There are many interesting questions concerning the applicability of usual 
physical assumptions and consequent results in the period between z ~ 10 15 
(kT m 400 GeV) and z w 10 10 (kT m 4 MeV). They are worth a brief paren- 
thesis, as they constitute the background to the discussion of the particular 
questions we shall be concerned with. 

First, concerning Particle Physics, it may be that the usual treatment 
of cross-sections need revision, due to scarcity of space. In that treatment, 
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particles are supposed to start in an interaction-free asymptotic region, in- 
teract in some finite intermediate region, and finish the process as a free 
object in yet another asymptotic domain. The ingoing and outgoing flows 
are compared to define the cross-section. Total volume is supposed to be 
much larger than the volumes of the particles concerned. Even if for decay 
rates and production processes with two initial particles (such as pair pro- 
duction) the volume factor cancels out ||, such processes can be inhibited by 
the smallness of complete (configuration plus momentum) phase space and 
by final state interactions. The usual elementary definition of the scattering 
matrix |9j involves a Dirac delta function in four momentum which actually 
assumes a very large volume. This is a rather difficult question, but one 
which can in principle be theoretically solved. 

Other, deeper Particle Physics aspects are to be considered. Nucleons 
are non-relativistic at z ~ 10 10 , but highly relativistic at z ~ 10 15 . How 
far will a proton, when energy grows from 4 MeV to 400 GeV, remain a 
proton ? Is it true that at the energies involved in the period we can already 
talk of deconfined quarks and unshielded gluons ? It is usually supposed 
that at high enough energies nucleons loose their identities and the system 



must be considered as a quark-gluon plasma fllQ"|. The energy at which that 



happens, however, remains unknown. The search for a signal of quark de- 
confinement, which has been actively looked for in the interval 15 — 200 GeV 
per proton, has not yet given a definitive answer [|ll|]. The experimental re- 
sults are consistent with the presence of deconfinement, but do not exclude 
other interpretations. This means that the signal is not unambiguous, is not 



conclusive ||12|| . What happens at still higher energies is simply not known, 



as the mechanisms of confinement and eventual deconfinement are as yet 



unclear even from the purely theoretical point of view | 13|| . We shall here 
suppose that protons keep their identities and examine the consequences. We 
shall find, as we proceed backwards in time, that radiation protons produce 
a pressure blockage at energies of the order of hundreds of MeV, much too 
low for confinement to take place. 

Concerning Thermodynamics and Statistical Mechanics aspects, the sys- 



tem should be treated, quite probably, as a finite system JT4|. Another point 



concerns fermions in general. The Pauli principle can be seen as a conse- 
quence of an effective repulsive potential between kin fermions, whose range 
is the thermal wavelength. This wavelength decreases with temperature. At 
the period under consideration, it will be very small for electrons, but large 
for non-relativistic protons. What has been said of phase space for cross 
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sections should be repeated in this context: pair production of protons, for 
example, would be inhibited by the presence of many protons in the final 
state. 

Anyhow, the first problem to be faced is much simpler. The picture 
we have of the early Universe comes from inserting ideal fluid equations of 
state in the Friedmann equations. In the PNS period, excepting the reactions 
leading to nucleosynthesis and their nuclei-breaking inverses, interactions are 
not taken into account at all, let alone the possibility of any abrupt behavior 
related to phase-transitions. Thus, the first thing to be done would be to 
consider a real, interacting gas. This is already a difficult enough task. We 
shall make a first attempt by taking into account the size of the nucleons. 
Protons and neutrons are, of course, the hardest known objects. Instead of 
structureless, pointlike particles, we shall assume a hard-sphere gas. 

A negative point is that the only known means to do it is by considering a 
hard-sphere potential, and potentials loose progressively their meaning when 
the particles represented become more and more relativistic. Taking into 
account pair creation and annihilation become more and more necessary to 
avoid difficulties akin to the Klein paradox ||. We shall do it as carefully 
as possible, in the hope of finding effects in the non-relativistic regime. The 
nucleon size-effect will be simulated by a hard-sphere with the nucleon radius. 
What we should take for the proton radius is as yet a matter of controversy 



15| , but the uncertainty between 0.80 fermi and 0.86 fermi is, of course, 
irrelevant for the gross estimate we have in view. We shall use 0.8 x 10~ 13 cm 
for the proton radius and undertake the usual backward path, from z ~ 10 10 
to z « 10 15 . 



4 A solvable model for primeval close packing 

We have presented in the introduction a very rough estimate of the red-shift 
at which a blockage can take place, and proceed now to a somewhat more 
elaborate evaluation. The result will be essentially the same as long as only 
the remnant protons are concerned. This more refined approach provides, 
however, general formulae which can be applied also to the protons belonging 
to the radiation. 

The hard-sphere gas is one of the great unsolved problems of Theoretical 
Physics. The virial coefficients have been calculated analytically only up to 
the fourth order (by Boltzmann). Numerical results exist for higher orders, 
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but the virial series has, if any, a poor convergence. The best picture of the 
system is given by a computer simulation refined through the use of Pade 



approximants. The outcome is a curve for the equation of state [jig] . It shows 
a clear phase transition (possibly two), in which the pressure grows steeply 
to oo. The curve can be parametrized by an equation of state of the type 

rih kT 
1 - n b /n c 

where n c = -^f, D being the sphere diameter. This equation can be quali- 
tatively understood in the "excluded volume" approach to the hard sphere 
gas. The canonical partition function for an iV-particle gas with interactions 
given by a potential Vy = V(|r$ — r^j) is 



1 f n\^ N ^ 1 \^ N 
Qn = J d 3 r 1 d 3 r 2 . . . d 3 r N e~' ^= 1 ^j=*+i 

. N-l N 

J d 3 n d 3 r 2 ...d 3 r N J] II e-^. 



X 3N 



=1 j=i+l 

For a hard-sphere potential, the integrand vanishes every time it happens 
that |rj — rj\ < D/2 for any pair of particles. Thus, it all amounts to 
forbid the regions with |rj— rj\ < D/2 for every pair, or to exclude the interior 
of all the spheres. As the potential simply vanishes outside the spheres, the 
final picture is that of an ideal gas in a volume reduced by the total volume 
of the N spheres. With the volume v s = | D 3 for each sphere, the equation 
would be p(V — Nv s ) = NkT, or p = nkT/(l — nv s ). There are actually 
some geometric factors before v s , as the excluded volume increases at close 
packing. In a configuration like 88 each sphere actually excludes a cube of 
size D. These geometric factors have been found by the authors of |16||, the 



second suggested phase transition corresponding to the still tighter packing 
of type & . 

Let us now apply to our problem. The proton radius gives D 3 = 
4.1 x 1CT 39 cm 3 , so that n c = 0.345 x 10 39 cm~ 3 . Using fllPP , we find 

— = 3.18x 10- u (1 + z) 3 n m h 2 . (59) 

The denominator in (|58|) vanishes when 1 + z — 3 i/3 10 2 1 / 4 3 , corresponding to 

fi b0 h 

kT « Jf 6 GeV . 
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Thus, ii h = 0.7 and Q b0 = 0.03, then kT « 300 GeV. The recently more 
favored value S7^ — 0.2 would give kT « 160 GeV^. This is an indication 
of close-packing for the protons existing today, but at the very high value 
z ~ 10 15 . The close-packing effect will actually take place at much lower z 
because of the pair-production effect, to be discussed in the next section. 
Notice anyhow that, even neglecting the pair-produced protons, the proton 
mean free path would be ~ (n c Ap) _1 , comparable to Xp itself. 
Using ( [40]), the baryon pressure ( |j)8|) will be given by 

9.2 x 10- 14 (1 + zf Q b0 T - — —JL ■ - - ^ . (60) 



c 2 p crit ' v ' ou 1 -3.18 x 10~ 44 (1 + zf n b0 h 2 

As there is thermal equilibrium at the period of interest, T will be the radi- 
ation temperature and consequently 

Pb n AO „ ln _i 3 o /i , .a 4 1 



2.48 x 10~ id Q 6 o (1 + z) 



c 2 p crit ' ' uv ' l-3.18x 10- 44 (1 + ^) 3 Q bQ h 2 ' 

The relation between this and the radiation pressure is 

^ = 3.17x 10- s Q b0 h 2 ' 



p 7 ' u l-3.18x 10- 44 (1 + z) 3 Q w h 2 

We can introduce the notation 

M = 2.48 x 10~ 13 n b0 h 2 , Q = 3.15 x 10~ 15 fi^ 3 h 2/3 (61) 

and write 

'"' M (1 + ^) 4 - ^ £ , . (62) 



c 2 p crit v y i-[g(i + z)] 3 

This is to be added to the Friedmann equation in the form (|55|), which 
becomes 

I | = 3/ + n^ + *^ + ^ T -^ J -A. (63) 



The solution is 



fc o A 



a o 



k A 
1 + 3 7 ° 



M 



3 J. (l-Q)V-i+Qx + QV V3Qx V3Q\ 

x < In — — == — h arctan arctan — > . (64) 

\ (l-Qx)Vl + Q + Q 2 2 + Qx 2 + Qj 1 ; 
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This is just fl5E|) , with the additional term proportional to M. As Q <C 1, a 
good approximation is 



/0) = fi 70 a; 4 - — a; 2 + - + x 3 



fc A 

1 + ^7- 77 



M 3 J yi + Qs + Q^ V3_Qx\ 

H x < in h arctan — — > . (bo) 

3 \ l-Qx 2 + Qxj K ' 

When x tends to Q~ x from smaller values, the functions f(x) and H(z) 
become infinite. The implicit solution for a(t) is still given by ([57]), with 
f(x) now given by flb~4]). 

We arrive thus at the following provisional picture. If we proceed back- 
wardly in time, there will be a value of the increasing red-shift for which 
matter pressure produced by the remnant protons becomes practically in- 
finite. This value corresponds to energies of a few hundreds of GeV. Of 
course, at such energies the very notion of potential is unacceptable. But 
this is only a first step in our reasoning chain. Actually, there will be a much 
larger number of protons in the medium. We have up to now neglected those 
which are produced in pairs by the radiation background. If that number is 
high enough, the pressure may become exceedingly high at lower energies, in 
which potentials still do make sense. 



5 Radiation-belonging nucleons 

The radiation background contains a large amount of massive particles as 
soon as the pair production 77 — > e~ e + threshold is attained. As to the 
process of hadron production 77 — > pp, it is clearly at work at kT around 
1 GeV. Actually, the reaction threshold is much lower [17], because of the 
huge number of photons. Even at lower energies, there are many photons 
with energy high enough to produce pairs. The number of radiation-created 
hadrons becomes of the same order of the number of photons a little above 
the threshold []nj , much larger in effect than the number of remnant protons 
considered previously. Let us say a few more words on these statements. 

The basic question would be: when we proceed backwards in time, from 
which energy on can we consider that the reaction 77 — > pp is in equi- 
librium ? This energy is important because, above it, the remnant protons 
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are negligible and the previous argument should be applied, instead, to the 
radiation-belonging protons. 

With an annihilation cross section o a ~ 2 x lCT 26 cm 2 , the annihilation 
mean free path of a proton is A a ~ — ~ 10 23 (1 + z)~ 3 cm. This gives ~ 

10~ 22 cm at 400 GeV, w 2 x 10~ 21 cm at 20GeV, and 10" 7 cm at 4 MeV. 
For the inverse reaction, pair creation, we could use the Breit-Wheeler cross- 



section [15 1, but the exact value is not necessary in our simplified approach. 
An estimate, using only mean free paths, can be made along the following 
lines. The mean free path for a photon, due to pair creation, is A 77 ^ P p = 
[n 7 o" 77 ^ P p] _1 . This means that, on the average, a 7 will meet another 7 
to produce a pair every time the volume A 77 ^ pp - o" 77 ^ P p is spanned by a 7. 
Thus, a traveling photon will "deposit" one antiproton at each volume of that 
value. But that volume is just l/n 7 , so that the concentration of antiprotons 
is roughly the same as that of the photons. 

A more precise description of the interplay between annihilation and pair 
production requires a detailed analysis of the kinetics involved. A kinetic 
estimate [Pj| gives r pa 1/44 pa 0.02 (corresponding to ~ 20 MeV) for the 



temperature above which there is equilibrium. We insist that, as soon as 
chemical equilibrium is attained, at kT pa 20 MeV, the number of antiprotons 
becomes enormous, of the order of n 7 . 

We can also estimate the temperature at which the number of pair- 
produced protons (or antiprotons) becomes larger than the number of rem- 
nant protons, by equating fl55| ) and (|5T|). The result of a numerical analysis 
is that the number of antiprotons overcomes that of the remnant protons 
at t pa 0.058. Thus, we can choose for security a reasonable value ~ 0.06 
(corresponding to ~ 60 MeV) and say that, for r above it, the reaction is 
in equilibrium and there will be a large amount of protons and antiprotons. 
At kT pa 60 MeV , these protons and antiprotons are nonrelativistic and for 
them the potential-based arguments are valid. 

We have repeatedly said that the concentration of protons is actually of 
the same order of magnitude of that of photons, many orders of magnitude 
above the number of the remnant protons we have considered in the previous 
section. In fact, using the numerical factor n 1 /rip = 4/3 discussed at the end 
of Appendix A, and equations ( |50| ) and (CT), we find 



^ = 2.9 x 10 7 [^h 2 ]- 1 . (66) 
rib 

There is consequently an abrupt jump in the concentration, and ([5TJ) will 



20 



change dramatically. To take this effect into account, it is sufficient to change 
the numerical parameters fl61|): M and Q 3 must be multiplied by given 



by (jog). Thus, the formulae of the previous section can be used, with the 
parameters M and Q of (^) replaced by 



M' = — M = 7.2 x 10- 6 ; Q' = ( ^ ) ' Q = 9.56 x 10~ 13 . (67) 

In that case, -fi— = M' (1 + -z) 4 i—f?v7TT m , or > as M ' = 0.65/i 2 ^ , 
_ 1.95/z 2 p 7 _ 0.65/?, 2 e 7 

C 2 Pc„ t ~ 1 " \Q' (1 + ^)] 3 C 2 Pcr, t ~ 1 " [Q' (1 + ^)] 3 C 2 p m / <(KNi 

The term VL l0 x in (|63|) is the radiation contribution. That equation becomes 

0.65 h 2 

f 1 - [Q> (1 + 

(69) 

The last term corresponds to the radiation contribution, with the term in h 2 
giving the nucleon interaction correction to the equation of state (0). The 
solution is still (|64]) , but with the replacements M — > M' and Q — ► Q'- 

Matter pressure now becomes infinite at z = 1.0 x 10 12 , or fcT = 232MeV. 
At such energies, it is possible that deconfinement, conjectured to happen 
at w 150Me^, has taken place. We are, however, neglecting some effects 
(discussed below) which would tend to lower that value. And even at ~ 
230MeV the non-relativistic argument, based on the notion of a potential, 
is valid. Notice that the energy at which pair production attains equilibrium 
is much lower: r w 0.05 corresponds to kT 1 Pd 50 MeV. A last point is that 
antiprotons attain a concentration comparable to that of photons a bit above 
the reaction threshold, which can be lower than the equilibrium temperature. 
We are, thus, quite probably overestimating the above energies. 



1/3 



:i+2) 5irfi) = 3/+ | (1+ ^- A+ ^ (1+2)1 



6 Conclusions and speculations 

The conclusion is that, as soon as pair creation starts up, the number of 
protons becomes so high as to create a blockage. This might mean that the 
reaction is inhibited by a "wall" of occupied phase space in the final state, 
or that the use of cross-sections is inadequate, or still that usual thermo- 
dynamics does not apply. It is common knowledge that, in consequence of 
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fundamental requirements, pressure in a fluid cannot go to infinity. It can 
at most attain the incompressibility limit, at which the equation of state is 
p = e and beyond which causality would be violated. With h — 1, this point 
would correspond to z ~ 10 12 . 

In its details, the primeval blockage we have dealt on has, unfortunately, 
a crucial dependence on the values of the as yet undetermined cosmological 
parameters. Some of the neglected aspects would, however, add to the re- 
pulsion and, consequently, to the effect. Notice that any effect leading to a 
higher effective value for the sphere volumes in (|58|) will bring the critical 
redshift to lower values. For instance, scattering theory will tell us that a 
sphere appears larger in the quantum case: naive shadow scattering by a 
disk of radius r will see a transverse area 4irr 2 , instead of the classical nr 2 . 
This would mean an effective radius twice as large, and a red-shift half that 
found above (z = 0.5 x 10 12 , or kT = 116 MeV). Still another possiblility 
is the breaking of Debye electromagnetic screening, which would lead to an 
increase of the electrostatic Coulomb repulsion. Finally, we have completely 
neglected the Fermi repulsion due to Pauli exclusion. This effect is feeble 
in the ultrarelativistic regime, but the thermal wavelengths ( f2~3|) or can 
be larger than the diameter of the excluded sphere for the range of energies 
found. 

We have not attempted to relate our results to the current negative- 
pressure approach to inflation pO, pTJ . We shall only indulge in a wild re- 



mark. As Lee and Yang have taught us long ago ]F2 , equations of state have 
the same analytical form both sides of a phase transition. If equation ( 58 ) 
holds for nb/n c > 1, the pressure will be negative "the other side" of the 
thermodynamic singularity. In effect, the expression for the total radiation 
pressure, including the nucleon-antinucleon pairs, 



p 7 



C Pcrit 



f2 7 o(l + z) 



1 0.65 h 2 

3 + 1-[Q' (l + z)} 3 



(70) 



or (using e 7 for the ideal energy density) 



Pi 



1 

3 + T 



0.65 h 2 



[Q'(l + z) 



+ 



0.65 h 2 



1 - 


"fcT 7 (MeV)~ 


3 




232 





-7 > 



will be negative in a short interval before the singular red-shift. A comparison 
of §L^) and (|52D shows that, for the large values of z we are considering now, 
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the energy density is dominated by the quartic term. Consequently, the ideal 
e 7 of the last expressions is the true radiation energy and (|7TD is indeed an 
equation of state. The equation would be of the exponential inflationary 
type p = — e type [^3| only for a very particular value of z. We can 
however, speculate on the possibility of a more general type of inflation. For 
a barotropic equation p = (7 — l)e an extended, power-law type inflation [25| 
occurs for 7 in the range < 7 < 2/3 ||24)| . For h = 1, this would correspond 
to a tiny interval 1.14 < Q' (1 + z) < 1.25 before the pressure singularity. 
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A Relativistic Gases 



We justify here some statements and formulae of the text, and rewrite some 
of the most usual expressions in units specially adequate to our case. 

To see how ( p3|) comes out, let us recall that the grand-canonical parti- 
tion function for a gas of non-interacting quantum particles with chemical 
potential /1 is the trace of the density operator: 



l(V,/3,ii)=tr 



E ( n onin 2 



V] V] e -/?(eo-M)no e -/3(ei-M)ni e -/3(e 2 -M)n 2 



IIE 



-/3(ei-/i)n 



n ni n 2 



IIE 



-P{t-n)n 



In this non-interacting case, each level contributes an independent factor. 
The system can have also internal degrees of freedom, which will likewise 
contribute separately. Suppose a single degree of freedom (spin, for example) 
taking g possible values. The partition function will be 



H(V,/3,//) = n 



n 



(72) 
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The kind of statistics appears in the summation, which is over the possible 
values of the occupation number n, from n = up to the maximum number of 
particles allowed in each state: 1 for fermions, oo for bosons. To treat bosons 
and fermions at the same time, we adopt the usual convention: upper signs for 
bosons, lower signs for fermions. Transforming the product into a summation 
by using the formal identity T7 e {- • •} = nJ ex P(l n {- • •})] = expJ2 e m {- • }j the 
above expressions lead to 

In ~ B > F (V, (3, M ) = «E In [l T e~^] . (73) 

It is convenient to use the fugacity variable, either the usual non-relativistic 
fugacity z = e^ 4 or its relativistic version Z = = ze^ mc . If we do not 
care about zero-energy states, the sum over the energy levels can be replaced 
by an integral over the momenta through the prescription J2e ~ * / d 3 xd 3 p, 
which leads to 

n 6 Jo L J 

Expanding the logarithm and collecting like terms, the partition function 
acquires the form 



E B > F (V, 0, z) = exp | ^ £ J d 3 pe -M( P ^ +m ^- m cZ] 1 (?4) 



gV ^ (±l) i_1 
i=i i 

The relativistic thermal wavelength (1231) appears now in the form 



¥J dpe ~Wfl' 

and the final expression for the grand-canonical partition function for a gas 
of non-interacting quantum particles is 

E-(V^He*p{^g<^^}, (76) 

or its equivalent 

-^(V, A ,) = exp L^^ft { ^fPK M) } . (77, 
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Here K 2 (x) is the modified Bessel function of second order. Limits can be 
found by using the properties 

jr l( ,)*yr e -( 1 + ^ + ...) forx »l; 

i^ 2 (^) ~ 2aT 2 ; i^i(x) « a; -1 for x << 1. 

Ki(/3mc 2 ) will appear only in the energy expression. The non-relativistic and 
the ultra-relativistic limits give ( p4|) and (|25|). The pressure and the particle 
number follow by standard thermodynamic relations: 

oo (_L.\l-l 1 

p V=H-h,a = 9 ATgLJ_ 1 ._ f _, (78) 



«E 1 



V,/3 e + 



_ ^ r tZ 3 p _ nF fv | y-i J 1 r7q ^ 

" /i 3 7 z -i e /VP*W?± 2 fe l J A£(Z/3) " 1 J 

The expressions in terms of integrals or of series are more or less convenient, 
depending on the application in view. We can extract the density number of 
particles at energy e, n e = g \z~ x e^ € =F l] -1 - The average energy, including 
the masses, is 

= 3pV + 4tt g ( ( V ^^z^e^K^ipmc 2 ) . (80) 
\ X c J V mc J j5 1 

The degree of degeneracy is 

^^~ 9 h { ) WW) ~ 9 h~^ e-^ 2 K 2 (Pmc 2 ) ■ (81) 

We are particularly interested in n = y. For a massless particle, the change 
of variables x = pc/kT can be used directly to give 

g t 3 f°° x 2 dx 
v 3 " 



n = ——t / 7 . (82) 
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For a gas of photons (using g = 2, z = 1), there are many expressions of 
interest : 



7 

n 7 — 



iV 7 2 r d 3 p 1 r 3 x 2 dx n ™ 1 



V h 3 J e^-l n 2 X 3 c Jo e»-l 2 §Af^(Z/3) ^ 

r 3 ~ 1 2 C(3) /A;T 7 \ 3 r 3 , s 

= 2^gp=^(^j = 0.244 ^. (84) 

The pressure is found to be p 7 = |® n 7 fcT 7 = E 7 /3V = e 7 /3. For a gas of 
fermions with g = 2 (like protons or antiprotons), 



AT 00 1 



3 




In the ultrarelativistic regime this becomes 

1 r 3 f 00 x 2 cfe 3 C(3) r 3 r 3 

n = = _ = 0.183 . 86 

p vr 2 X 3 C Jo e x + 1 2 7T 2 A^ A 3 . v ; 



A factor ^ = | comes from the fermion repulsion effect, encapsulated in 
the sign in the integrand denominator, opposite to that in (|83|). It is more 
difficult to pack fermions than bosons together. This can be seen also from 
the limits of the degeneracy index flST]). In this ultrarelativistic regime, its 
values are, for photons and antiprotons, respectively, n 7 Afj R = 0.244 n = 
0.766 and npuR A.fj R = 0.183 tc = 0.575. The relativistic or nonrelativistic 
character of the protons depend on the above numerical factors. They are 
possibly irrelevant to rough estimates of the main text, but may come to be 
important in a more detailed consideration of Fermi repulsion. 
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